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I ■ The information content and properties of the cross section for atom scat- 

tering from a defect on a fiat surface are investigated. Using the Sudden 

& : 

approximation, a simple expression is obtained that relates the cross section 
to the underlying atom/defect interaction potential. An approximate inver- 
O ■ sion formula is given, that determines the shape function of the defect from 

the scattering data. Another inversion formula approximately determines the 
potential due to a weak corrugation in the case of substitutional disorder. 
An Optical Theorem, derived in the framework of the Sudden approximation, 
plays a central role in deriving the equations that conveniently relate the in- 
teraction potential to the cross section. Also essential for the result is the 
equivalence of the operational definition for the cross section for scattering by 
a defect, given by Poelsema and Comsa, and the formal definition from quan- 
tum scattering theory. This equivalence is established here. The inversion 
result is applied to determine the shape function of an Ag atom on Pt(lll) 
from scattering data. 
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I. INTRODUCTION 



The pioneering work by Poelsema, Comsa and coworkers |]T]-§| on atom-surface scatter- 
ing in the last decade, has led to the introduction and application of the concept of the 
cross section of a defect on a surface in studies of questions related to surface disorder. 
Extensive experimental work in the field (as summarized in ref. 0) has focused primarily 
on the measurement of the dependence of the cross section on the surface parameters, usu- 
ally the coverage dependence. Some efforts have also been directed to the dependence on 
the scattering parameters, such as incidence energy and angle. These cross section studies 
have attracted significant theoretical interest, ]|7|-|l4j] concentrating almost exclusively on the 
dependence on the scattering parameters. Relatively little attention has been paid to the 
inverse question of the information content of the cross section. Most work in this direction 
has been experimental, again focusing on the use of cross section data to study coverage- 



dependent questions, such as the onset of Ostwald ripening, [15] and whether growth proceeds 



in a layer-by-layer fashion. pq-pT| With very few exceptions, [|14] , p2|| most studies to date 
do not deal with the issue of the structural information contained in the cross section. For 
instance, can the cross section be used to extract the shape of an object on the surface, or 
its electron density profile? Can the interaction potential between an atom and a surface 
defect be obtained from the cross section data? In the present paper, we address these 
questions theoretically, by focusing on what can be learned from studying the cross section 
of an individual defect on a surface. 

To make contact with experiment, through the operational definition of the cross section 
proposed by Poelsema and Comsa, || it is required to show the equivalence of the latter 
definition to the standard formal one from scattering theory, which we shall employ. Here 
we will establish this equivalence. An analytical study of the cross section is most easily 
performed by relating the cross section through an optical theorem to the scattering am- 
plitude. This has been done by several authors. Hl8|,|7|,|9| We derive an optical theorem for 
atom-surface scattering by use of the Sudden Approximation (SA), |23| which has two sig- 
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nificant advantages: it leads to a simple explicit form for the cross section, and is expected 
to be fairly accurate, as the SA performs best for the calculation of specular amplitudes. By 
employing our optical theorem result, we will answer a number of specific questions related 
to the information that can be extracted from the single defect cross section. 
First we shall examine the relative role of the long range dispersion forces vs the short range 
repulsive part of the potential in determining the interaction of a He atom with an adsorbed 
defect. It is a well-recognized fact in surface scattering theory, [|4],[?],|ll],|! that the magnitude 
of the total cross section of a defect is dominated by small angle deflections caused by the 
long range forces. We will show that they are mostly effective in setting the magnitude of 
the cross section, whereas at least for the high-energy regime the short-range forces tend to 
determine its qualitative shape. 

A second issue that will be pursued is whether the cross section be used to perform an 
inversion of the He-defect interaction potential. An affirmative answer will be given in two 
important cases: the cross section can be used to invert the semiclassical shape function of 
a defect on a surface, and to approximately invert the interaction potential of a He atom 
with a weak corrugation, e.g., due to dilute substitutional disorder. The results are applied 
to obtain the shape function of an Ag atom on Pt(lll), using cross section data from P. 
Zeppenfeld. f25[ 

The structure of the article is as follows: we discuss the standard and operational definitions 
of the cross section in section [H], and prove their equivalence. In section [TTJ, we briefly review 
the SA, and derive the optical theorem for atom-surface scattering within the SA. We then 



proceed to discuss a number of applications of the optical theorem, in section [IV| The first is 
a discussion of the influence of the short range repulsive part of the atom-surface interaction 
potential. The second is perhaps the central result of this paper: an inversion scheme for 
the shape function of an adsorbate on a surface from cross section measurements. We apply 
this result to data for Ag/Pt(lll). The third application is the approximate inversion of 
the potential from the cross section data. Concluding remarks are presented in section [V] 
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II. DEFINITION OF THE CROSS SECTION 



Consider a scattering experiment of an atom beam with wave- vector k = (k x , k y , k z ), 
impinging on a flat surface (the z = plane). Let R = (x, y) denote a vector in the plane. A 
typical experimental setup will involve a highly monochromatic beam with a lateral spread 
in the surface plane. This beam is suitably described by a set of states \^i n ) = |$r). 
After impinging upon the surface, the beam is scattered, in part by the defect whose cross 
section we are interested in. To derive the cross section, we start from the definition of 
the quantum differential scattering cross section (QDCS) for two-body collisions. |2B| This 
definition generalizes quite naturally for atom-surface scattering, as we will demonstrate. 
Let w(dQ <— $r) denote the probability that a particle incident upon a target as a wave 
packet $, displaced laterally by R with respect to the origin, emerges after an elastic collision 
in the solid angle dQ. Then the average number of observed scatterings into dQ is: 

N sc (dn) = J>(dfi <- $hJ « J d 2 Rn inc w(dtt <- $r) 

i 

where n inc is the incident density, which we assume to be uniform over an area A, i.e., 
n inc = N inc /A, so that: 

N sc (dn) = (1) 

where da/dQ is the differential cross section, defined as: 



da 
~d£L 



dn = a{dVL)= f d 2 Rw(dtt+-$ R ) (2) 

J A 



Thus, as is well known, the QDCS is seen to be a measure of the target area effective in scat- 
tering the incident beam. It must be emphasized that the definition of the QDCS excludes 
the forward direction. Classically speaking, the reason is that one cannot distinguish a for- 
wardly scattered particle from one that did not interact with the target at all. The analogue 
in surface scattering from a defect, is that one cannot distinguish a specularly scattered par- 
ticle from one that was scattered by the flat part of the surface. More formally, the reason 
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is that inclusion of the forward direction would lead to a divergence of the cross section. 
We would now like to discuss the connection between the latter definition, which is widely 
accepted for two-body processes, and the operational definition proposed by Poelsema and 
Comsa for atom-surface scattering, and applied extensively for cross section calculations 
from diffraction data. Their formula reads: 

S op = -lim^^ (3) 
n^o 1 an 

where S op is the total scattering cross section, I s , I represent the specular scattering intensity 
from the corrugated and flat (smooth, defect-free) surfaces respectively, and n is the number 
of defects per total surface area. Xu et al. |10| employed this formula to derive an expression 
for the cross section based on Gaussian wave packets, and showed that their expression can 
be interpreted as an optical theorem. We will show that the operational definition (eq.([3])), 
is in fact identical to the standard definition of quantum scattering theory (eq.(^)), without 
recourse to a specific method of calculating the quantities that appear in the formulas. This 
is true, provided the proper conditions are established to relate two-body and atom-surface 
collisions. The conditions under which the QDCS is defined, are: 



• the incident particles have a sharply defined momentum po and are randomly displaced 
in a plane perpendicular to p , 

• the distribution of scatterers must be such as to avoid coherent scattering off two or 
more centers. 

When interpreted in terms of atom-surface scattering, we see that the first condition 
may be traded by random lateral displacements parallel to the surface plane (maintaining 
a sharply peaked momentum). The second requires that the defects are distributed very 
dilutely. This is, however, exactly the assumption underlying the operational definition, 
and of course holds for scattering off a single defect. Thus the conditions assumed for 
the definition of the two-body cross section are naturally extendible to the atom-surface 
case, and suit the operational definition. In order to demonstrate the equivalence of the 



definitions, we first replace N SC (Q)/Ni nc by I(dQ)/Io dQ in eq. ([!]), where / is the q/j^-specular 
intensity. Iq is the incident intensity, which is clearly equal to the total scattered intensity in 
the case of a flat surface, and is hence indeed the same quantity as defined in the discussion 
of the operational definition above. With these replacements, the total cross section is: 

E = I ^ dn = T I J ( fi ) dn ( 4 ) 
J dQ io J 

Next consider taking the limit in the operational definition. This can be done by assuming 
that only one defect is present on the surface, so that for all practical purposes dn = 1/A, 
A being the total surface area, and dl = I s — I . Then from eq.(|[): 

S op = ^(/ -/ S ) (5) 
-»o 

But clearly, I — I s = J I(Q) dQ, so that S op = E, as claimed. 

III. OPTICAL THEOREM FOR SURFACE SCATTERING IN THE SUDDEN 

APPROXIMATION 



We very briefly review the SA. For more details see e.g. the review by Gerber. [23 



A. The Sudden Approximation 

Basically, the SA requires that the momentum transfer in parallel to the surface be small 
compared with the momentum transfer normal to the surface, i.e. |p7| , p8| 

|q' - q| < 2k z (6) 

where k z is the incident wave number in the z direction, q is the incident wave vector in 
parallel to the surface plane, and q' is any intermediate or final wave vector in parallel to the 
surface plane which plays a significant role in the scattering process. The S-matrix element 
for scattering from q q' may then be shown to be given by the function: 

(qlSK) = T / dRe 2 ^ (R) e iAq ' R (7) 

A J surf 
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where R = (x,y) and f](R) is the scattering phase shift computed for fixed R, given in the 
WKB approximation by: 



roo 

nCR) = / dz 

Jz t (R) 



2 



k 2 2mV(R,z) 



k z z t (R) (8) 

'zt(R) 

Here Zt(R) is the classical turning point for the integrand in eq.(|8|). The normalization 
is chosen such that the unit-cell A is the average area available for one adsorbate, i.e, is 
determined by the coverage. This choice is suitable for low coverages, when the simultaneous 
interaction of a He atom with two neighboring adsorbates is small. Of course, an adsorbate 
may "spill over" beyond its unit-cell due to the long range of the potential, and hence the 
integration in eq.(|7]) is over the entire surface. Condition @ for the validity of the Sudden 
approximation is expected to break down for systems of high corrugation. For instance, 
certain isolated adsorbates on an otherwise flat surface can represent, for realistic parameters, 
a very substantial local corrugation. Nevertheless, previous calculations have shown that the 
Sudden approximation reproduces rather well many features of the scattering from isolated 



adsorbates. |29| Features for which it breaks down are, e.g, intensity peaks due to double 
collision events, in which the incoming atom first hits the surface and then the adsorbate (or 
vice versa), which are a particularly sensitive manifestation of a strong corrugation (we note 
that a double-collision version of the Sudden approximation has recently been developed 
301). However, due to the assumption @, the Sudden approximation is particularly useful 



for the evaluation of specular intensities at incidence angles close to the specular direction, 
as confirmed by comparison with numerically exact wave packet calculations. Jll]| In the 
context of the present work, only the specular intensities are required (which are obtained 
by setting q' = q in eq.(|7])), which constitutes the most favorable condition for the Sudden. 



B. Derivation of the Optical Theorem 

The ensuing discussion closely follows that of Taylor [^6| for two-body scattering. We 
proceed to derive the optical theorem for surface scattering within the SA. 
The probabilities w(dQ <— 0r) in eq.(0) are: 
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roc 

W (dn^v in ) = dn P 2 d P \q out (p)\ 2 ( P = hk) (9) 

JO 

since w(d 3 p <— \I/ lTl ) = d 3 p \^ out (p)\ 2 is the probability that long after the collision the 
particle incident as emerging as ^ out , has momentum in (p, p + dp). Only the direction 
of p is of interest, as we are assuming elastic scattering. The "in" and "out" states are 
related by the S-matrix: \^ ut) = S \^i n )- In the momentum representation: 

*««t(p) = / d 3 p' {p\S\p') % n {p') (10) 

The combination of eqs. (PfflTOl) yields the differential cross section. Following the discussion 
in the preceding section, it can be seen that there is nothing which essentially restricts the 
use of eqs.(3,p|-|T0|) to two-body scattering, and with the restriction of scattering angles to the 
half-space above the surface (z > 0), they may as well be applied to atom-surface scattering. 
We now proceed to calculate the differential cross section within the Sudden approximation. 
Let us denote p = fo(q, k z ), where q is the wave- vector component in parallel to the surface, 
and Aq = q' — q. We require the familiar "scattering- amplitude" /, in terms of which the 
Optical Theorem will be expressed. To find its form within the SA, we first write down a 
general expression for the S-matrix element, in the case of scattering from a surface with a 
defect. This expression must be composed of two terms, reflecting specular (the analogue of 
forward scattering in the gas phase) and off-specular scattering. The first arises in the case 
of scattering from the flat part of the surface, whereas the second is due to scattering from 
the defect. Energy is assumed to be conserved, so the general expression for the S-matrix 
element (pl^p') contains an on-shell delta-function of the energy E p = p 2 /2m: 

(p\S\p') = 6(p z + p' z ) <y(ftAq) + ^^S(E P - E p ,)f(p «- p') (11) 

The identification of the scattering- amplitude by comparing eqs.(|7]),([ll]) is now made as 
follows. First, the phase-shift 1] is constant for scattering from a flat surface, and may well 
be chosen zero. In this case the SA yields (q| 5"|q') = £(Aq) (eq.(^) with r/(R) = 0), in 
agreement with the general form for specular scattering. Second, we consider the remainder 
after subtraction of the specular part from the RHS of eq.(0): 
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3( q <_ q ') = ^ f dR (e 2ir >™ - l) e iAq R (12) 

Z/l ./surf ^ ' 



so that: 



(q|5|q / ) = ^(Aq) + ^(q^q / ) (13) 
Energy conservation is implicitly assumed in the SA expression eq.(^). Still, comparing 



eq.([13D and eq.flllD, g is not quite the scattering amplitude yet, since it does not include 
the action of the S-matrix on the wave-function part in perpendicular to the surface. The 
SA derivation [^7| yields only the action on the part parallel to the surface. However, we 
can obtain the corresponding factor by comparison to the well known partial- wave expansion 
result. If we assume a plane wave He-beam at normal incidence to a perfectly flat surface with 
a spherically symmetric He-adsorbate interaction, the partial-wave expansion conditions 



apply. The gas-phase expression for the scattering amplitude is then: |26| 

/(k ^k') = ^E r r(k) \e 2i5 ^ - l] F™(k')* (14) 

where YJ m (k) is the spherical harmonic whose argument k denotes the polar angle (6,<ft) 
of k, and 5i(E p ) is the phase-shift in the angular momentum basis {\E,l,m)}. Comparing 
eqs.( p^) ,(|l4D, we see that the contribution of the component perpendicular to the surface is 
a factor of 1/k, so that the SA scattering amplitude is: 

/(p^p') = ^(qW) (is) 

We emphasize at this point our strategy for deriving the optical theorem: the analogue 
to the gas-phase forward vs non-forward scattering, is the flat-surface vs defect scattering. 
Within the SA, this distinction is easily made by separating the vanishing part from the 
non-zero part of the phase-shift. 

We are now ready to calculate the differential cross section within the SA. Substituting 
eq.flTTD into eq.(|T0D yields: 

*^(p) = <Mq, ~k z ) + / d 3 p' 5(E P - E p/ ) /(p <- p') % n (p') (16) 

9 



The first term is the wave scattered from the flat part of the surface, the second by the 
defect. Next we restrict observation to defect scattering only. Using the lateral displacement 
assumption from section || we then obtain: 

*«*(p) = ^ / dV 5(E P - E p/ ) f(p <- p') e^ R P'$( P ') (17) 

Substituting eq. ([T7|) into eq.@ and then into (@): 

2irm z fi J A Jo 
J d 3 p' 5(E P - E pl ) /(p «- p') e-i R p '$(p') x 

| dV (5(-E p - E pll ) (p <- p") e^P>(p") (18) 

The R integration yields (27r) 2 5(q // — q'). It is easily checked, using some well known 
properties of the delta-function, that consequently: 

J d 2 R e ^ R (p"-P') 5(E P - E pl ) 6(E P - E v „) = (2tt) 2 ^5(p" - p') 6(E P - E pl ) (19) 

Inserting this into eq . (|18D and applying the momentum delta-function, we obtain: 

hijV) roo r 1 1 

a{dQ) = ^hl p2dp /^V F ^-^0b(q-q0l 2 l*(p0l 2 7^i ( 20 ) 

Next we use S(E P — E p >) = —8{p — p') to do the radial integration over p: 

a(cM) =dnj d 3 p' _L \g(q <- q') $(p')| 2 (21) 

where |p| = |p'|. To simplify the last integral and obtain the standard form, we recall that 
$ represents the incident wave-function (at zero lateral displacement). The simplification is 
obtained by the physically plausible assumption that $ is sharply peaked about the incidence 
momentum po- More quantitatively, we assume that jj-yQ (q <— q') does not vary appreciably 
in the region of $'s peak, so that it can be taken outside of the integral sign, and be replaced 
by its values at k . Using further the normalization condition / d 3 p' |$(p')| 2 = 1, we finally 
obtain the familiar looking result for the differential cross section: 

a(dn <- ko) = dSl-±- |<?(q <- q )| 2 (22) 
10 



where from eg. (|T^) : 

g (q <- q ) = /" dR sin (77(H)) e ir?(R) e * (q °- q) - R (23) 

What we are after is the total cross section for scattering by a defect. We next employ the 
differential cross section result to express this in terms of an optical theorem for surface 
scattering. To do so, we use the manifest unitarity of the SA expression, eq.(^). Separating 
once again the specular and non-specular parts of the S-matrix, now in operator form: 
S = 1 + R, and applying the unitarity condition, we have: R + R' = —RR^. Inserting a 
complete set and taking matrix elements: 

(p'IaIp) + (p|iV>* = - / dh"{v"\R\v'Y <P"|fl|P> (24) 

But it is clear from eq.([TT|) that: 

(p'| J2|p) = — p(q <- q') 5{E p - E p ,) 
Inserting this into eq.(p4|) and factoring out a common 5: 
- (g(q[ «- q) - </(q - q')*) = ^~ [ d 3 p" 6(E p - E p „) gtf <- q')* <7(q" - q)-L (25) 

p 271171 J VP) 

The optical theorem results by considering the specular scattering (q = q') and integrating 
over p". In contrast to two-body scattering, in the atom-surface scattering case only the 
space above the surface is available, which can be taken care of by a factor of |. Then, using 
/ °° dp' h(p') S(E p — E p /) = ^h(p) and performing the radial integration: 

/ dY "' = \i ^ r° /2 dp " s{Ep - Ep,,) ^ (q,/ ^ q,) 9{q " ^ q) ^ 

Using this and the specularity condition in eq.(|2"5p, we obtain: 

*Im fo(q <- q)] = ±-±- J d^, \g{c[ <- q)| 2 (26) 

Comparing the RHS of eq.(^) with that of eq. (|2"2"D , we observe that it is just the total cross 
section for scattering at incidence momentum p, i.e.: 
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S(p) = J dQ p , a(dQ p , <- p) = J dQ p , \g(cf «- q) 
whence finally, the Optical Theorem: 



^k) = ^imfo( q <-q)] (27) 

From eq.(p3[) we obtain an explicit form for the total cross section: 

s(k) = Hi:I,/ Rsm2 " (R) (28) 

Equation ( ffijf) is the optical theorem for surface scattering within the SA. As a final remark, 
we did not consider the effect of Bragg scattering from a periodical underlying surface. As 



shown by Xu et al., [10] this contribution should be counted along with the specular ampli- 



tude, and its incorporation in the present framework should pose no particular difficulty 

IV. APPLICATIONS 

We now illustrate a number of applications of our result. 

A. Solvable Model of a Single Adsorbate 

It has been known for long that the long-range attractive interaction is dominant in 



determining the He- adsorbate cross section. [24,OjqJ3] We wish to identify the contribution 



of the short-range repulsive forces to the cross section. For this purpose, we will completely 
ignore the long-range attractive forces and consider a hard wall He-adsorbate interaction, 
where the adsorbate is described by a shape function £(R). Then the potential takes the 
form: 



V{r) 



: z > £(R) 
oo : z < £(R) 

From the expression for the phase- shift eq.@, it is clear that in this case: 

t7(R) = -k£(R) (29) 
12 



We now consider a fully solvable model of a single adsorbate, which although highly artifi- 
cial, will provide us with valuable insight regarding the role of the repulsive interaction in 
determining the general nature of the cross section. Let us assume a cylindrically symmetric 
shape function of the form: 

: a > 1 

e(R) = /»C(|R|/0 C(«)= (30) 

I 1 - a n : < a < 1 

where h is the height of the adsorbate above the surface (effectively the strength of the 
coupling between the He and the adsorbate). I = 7 L (0 < 7 < 1) is a characteristic range of 
the He-adsorbate interaction. The unit-cell area (see section p 1 111 A] ) is taken as A = 7rL 2 . 
The linear extent of the entire surface is denoted R s . We will demonstrate this choice of C(°0 
to be analytically solvable in terms of known functions for n — 1, 2, 4. For n = 1 one obtains 
a cone, for n = 2, 4 the shape is a convexly deformed cone. It has a physically unreasonable 
sharp edge for a = 1, yet the cases n = 2, 4 should be a crude, but reasonable model of the 
main features of the adsorbate shape function. In order to check the influence of the sharp 
edge, we may consider instead of ( a concave shape function, which trades the edge with a 
sharp tip at a = 0: 

: a > 1 

£(R) =hK{\R\/l) K{a) = \ (31) 

I (1 -a) n : < a < 1 

This model is solvable for n = 1,2. 

Using the symmetry, we have from the expression for the cross section within the SA 
(eq.([28|)), for the convex case: 

E(k) = , - 1 / da a sin 2 [h k z C(a)l (32) 
k k z Jo 

where e = R s /l. Defining (3 = hk z , and letting I^iP), In(fi) st an d for the integrals (without 
the prefactors) in the convex and concave cases respectively, they evaluate to: 

/<(/?) = = e - + C ° sm ~ C ° S(2/?(1 " e)) ■ esin (2/3(l-<0) 



4 8P 2 4 j3 
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c m _ e 2 sin(2/3) - sin (2/3(1 - e 2 )) 



4(/3) 



F c -2J/3/7T -F c [2J/3/*(e-l 



4 8/3 

sin(2 /3) - sin(2 (5{e - 1] 
V 8/3 



,2 



1 Fk 



4(/3) = - - - p (cos(2/3) F c [2e 2 ^/tt] + sin(2/3) F s (V ^/tt] J (33) 

where F c , F s are the Fresnel integrals: 

f x ( 7ry 2 \ f x ( ny 2 \ 

F c (x) = J o cos I — I dy F s (x) = sin I -y- I dy 

An analysis of the models in the simple case of e = 1 yields the following results (see figure 
[T] - cone, figures @, ^ - convex cases, figure [|- concave cases): 

• In both the concave and convex cases the cross section generally decreases with in- 
creasing incidence wave-number k z . 

• The convex model exhibits more noticeable oscillations. These are due to interference 
between He particles striking the top of the adsorbate and the surface. Since the top 
has a larger surface area in the convex case, the increased oscillations are expected. 

• The cross section is generally larger in the convex case. This is also expected due to 
the larger surface area of the top. 

• The cross section has a finite limit for k z — > 0, which can easily be found by a first- 
order Taylor expansion of the sin 2 in eq. fl3"2"| ) . This limit is as a matter of fact beyond 
the region of validity of the SA. However, relying on the impressive success of the 
SA under unfavorable conditions of significant corrugation, [2{| we consider the limit 



anyway: its dependence on the corrugation parameter h and the concavity /convexity 
parameter n is of high interest. We find: 

, /r , 32vr 2 h 2 n 2 

E c (k 0) = — (convex) 

1 ; (1 + n) (2 + n) V ; 

32tt 2 /i 2 

T, K (k — > 0) = — r (concave) 

V ; (l+7i)(l+2n) V 1 



14 



The corrugation parameter h is thus seen to be of major importance in determining 
the cross section, which in the limit of low incidence energy scales as its square. The 
dependence on the exponent n is, as expected, opposite in the concave vs. convex 
case (see figure |5|): in the former, the cross section monotonically decreases with n, 
whereas it increases in the latter. This is once again due to size of the surface area of 
the adsorbate top in each case. 

Due to the absence of the attractive part in our model, the resulting cross section values 
are too small. Typical values for e.g. CO/Pt(lll) |2||| are in the range 100 — 300 A 2 . 
This deficiency of the hard wall model was already observed by Jonsson et al. p4|,[7| 



• The integrals In(P) tend asymptotically to a constant value with k z . Therefore 

the decrease in the cross section is entirely (apart from oscillations) due to the k~ 2 
factor. 

In conclusion, it appears that the cross sections we obtained exhibit all the general features 
of the experimentally observed ones: the oscillations, the decrease with increasing incidence 
energy and the finite limit for vanishing incidence energy. We conclude that while the long- 
range attractive part is important in the setting the magnitude of the cross section values (the 
stronger the interaction, the larger the cross section), at least for the high-energy regime the 
energy dependent features are already determined by the short-range repulsive interaction. 
In the next section we will significantly simplify the relation eq . (j32[) between the cross section 
and the shape function, without any further approximations, in order to be able to eventually 
invert the shape function from cross section measurements. 

B. Inversion of Shape Function from Measurement of Cross Section 

In this section we will demonstrate how the optical theorem result may be used to 
approximately solve the inversion problem for the He-adsorbate potential. The SA was 
already employed by Gerber, Yinnon and coworkers [3"I,3"2|| for an approximate inversion of 
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the potential, by using the full angular intensity distribution. Here we will demonstrate a 
less general result for the inversion of the adsorbate shape function, which will, however, 
have the advantage of relying on a simpler measurement, of the specular intensities alone. 
We consider the case of a cylindrically symmetric adsorbate. If our unit-cell and the surface 
are circles of radii L and R s respectively, then from eq.(P5[) we have for the adsorbate cross 
section: 

32tt 2 „ r R ° . o „ 1 



Z(kz,0in) = —pn- 2 tt-^tt / sin 2 V (R)RdR = — -R 2 - / cos (2 V (R)) RdR 



(34) 



where 

I? cosf 



a 



32vr 2 

(We note in passing that the often discussed [0] cos(6'j n ) factor, in being the angle between 
the incident direction and the surface normal, arises naturally in our formalism). Our 
demonstration in section |I| that this cross section is equivalent to that measurable by a 
specular attenuation experiment, ensures that the last equation can be regarded as a simple 
integral equation for the phase-shift r\ in terms of measurable quantities R s , L and E(/c z , 0j n ). 
The measurable part is: 

M{k z ) = l -R s 2 - ak 2 z Z(k z ,9 in ) (35) 

Let us now make the hard wall assumption again, which leads to the form of eq. (|2~9D for 
the phase-shift. By doing so, we will be able to obtain the form of the shape function £(R) 
of the adsorbate. Of course the concept of a shape function is rather artificial, yet it is of 
great interest in characterizing surface roughness, and is not unrelated to what is obtained 
in STM measurements. Defining for convenience: 

e0R)=2£0R), (36) 

we can express the integral part as: 
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'6 



[ Ra RdR cos (k z Q(R)) 
Jo 



(37) 



so that eq. ([3"4"D becomes: 



Is = M(k z ) 



(38) 



We next assume 9(i?) to be a monotone, single- valued function, so that we can consider its 
inverse i?(0) (see figure Thus: 



e 



/e(o) <i6 
where we defined: 



1 /-©ORs) / 
i2(6) cos(A;,e) = -= / dS F(Q) (e ik * 6 + e 
V27T ie(o) v 



■ik z e 



F(6) 



27rrf (i? 2 (9)) 
~dQ 



R(Q) is confined, so it is convenient to define a "box" function: 



6(6) 



1 : Q(R S ) < 6 < 0(0) 
: else 



which allows us to write Iq as a Fourier transform of a product: 



Ie = F[F(Q)b(e); k z \ + ^ [F(6) 6(6); k 



where we introduced a notation for the Fourier transform, 



^ [/(*); y] 



V27T •/— oo 

To isolate F(Q) we apply an inverse Fourier transform to eq. (f42|) , yielding 
F(6)6(6) + F(-6)6(-6) = F' 1 [I e (k z ); 6] 



(39) 



(40) 



(41) 



(42) 



(43) 



But > so that by definition of the box function the second term in eq. ([43D vanishes, 
and we obtain from it, upon inserting the definitions of M (eq. fl35|) ) and F (eq. fliOD ): 

d(R 2 (Q)) 



dQ 



ttR s 2 5(0) - olT" 1 k 2 z Z(k z ,e in ); 6 for Q(R a ) < 6 < 9(0) (44) 



17 



If we assume that the shape function never fully vanishes (Q(R S ) > 0) then 5(0) = 0, so 
that finally: 

ffll^^fe ft); 91 (45) 

rfe 87r 2 \/27r L J 

The last equation is the sought after inversion. Its RHS should be considered as calculated 
from the experimental data on the cross section as a function of the incidence wavenumber, 
at a given incidence angle. The inversion is completed by solving for i? 2 (B) under the 
condition R(0) — > oo, and retrieving the shape function £(R) = |0(i?). 
In practice, however, it is more reasonable to assume a functional form for Q(R) with free 
parameters, to be fitted based on the experimental data. The somewhat artificial concept of 
a shape function may not justify an attempt to determine an actual function from the cross 
section data, as implied by eq. (|45|) . This could be problematic due to the discreteness and 
finiteness of the available data set. If we assume this approach, a further simplification can 
be obtained by returning to eq.(^) and integrating by parts (with Q(R S ) = - an excellent 
approximation) : 

f 0(O) 



R 2 s -k z dOR 2 (Q) sin(fc 2 6) 
Jo 



(46) 



Equating this according to eq . (|38|) to M(k z ), and assuming that we have guessed a functional 
form for Q(R) with a set of free parameters {a^}, the inversion result may now be stated in 
the form: 

( °' n -Y, ein (k z ) = g(k z ; {ai}) g (k z ; {c^}) = [ dB R 2 (Q; {aj) sin (k z 6) 

L z K z Jo 

(47) 



(4tt) 



The remaining task is to find the set {a{\ which gives the best fit for the chosen form of g to 
the cross section data. In the next section we shall implement this approach on experimental 
data. Alternatively, the last result may be viewed as a theoretical expression for the cross 
section, within the SA and a hard wall model. 
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C. Fit of Shape Function for Ag on Pt(lll) 



We now employ our inversion result to obtain a fit for a shape function of an Ag atom 
on a Pt(lll) surface, based on He scattering measurements by P. Zeppenfeld. [[25| The 
shape function is a legitimate and highly relevant object of study in the field of surface 
roughness. It basically corresponds to an equipotential surface, and therefore within a certain 
approximation, also to the core electron densities. As a result, it is a complementary 

quantity to what is measured in STM experiments, where one probes the electron densities 
at the Fermi level. Assuming cylindrical symmetry, the shape function should have the 
general form depicted in figure ^ it should be flat and smooth at its top, and decay to zero 
height far away from the nucleus. This leads us to consider two simple forms for the shape 
function: a Lorentzian and a Gaussian, with two free parameters I, 0o each. Thus: 

&l(R)= e ° x2 or G ( J R) = 6 o e-(T) 2 (48) 



1 



Inverting these for i? 2 (9) leads to integrals g (k z ; {«i}) (eq. (|j])) expressible in terms of 
known functions: 

n fe {»■ e,}) = £ ^ k *^ ~ ») + ggjjkgg)) (49) 
sg (*,; {/, e„» = ~ (r + in(fe e„) - a(k, e„)) (so) 

where si(x),ci(x) are the sine and cosine integrals respectively, and V is Euler's constant 
(0.5772..). Both and go contain an oscillatory part, expressing the interference between 
parts of the He beam striking the top of the adsorbate and the flat surface, and both decay 
as k~ 2 for large k z . The parameters I (the He-adsorbate interaction range) and L (the unit- 
cell extent) appear only in the combination l 2 /L 2 , and may hence be treated as a single 
dimensionless parameter 7 = l/L. This result is easily seen to hold for any shape function 
of the form O (R/l), as was already demonstrated in section [1V[1V A| The dependence of the 



cross section on the He-adsorbate interaction range thus contains no surprises: the cross 
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section scales as the square of this range. It is the adsorbate height Go, i.e., the coupling 
coefficient between the He and the adsorbate, which plays the interesting role. It determines 
the frequency of oscillations of the cross section, and affects its magnitude. 
Using the results for the Lorentzian and Gaussian models, eqs.(^,^) respectively, in the 
expression for the cross section, eq.(|4"T|) we calculated best-fits for the parameters 9 and 
7. The fits are shown along with the empirical cross section values |25] in figure The 
deviations are 1.47% for the Lorentzian, 1.14% for the Gaussian. The fits do not follow the 
oscillations in the data closely, but these are uncertain anyway due to the large experimental 
error (~ 20%). The values found are: 
Lorentzian: O = 1.5lA, 7 = 0.93 
Gaussian: 6 = 1.05A, 7 = 1.83 

The values for the adsorbate "height" 0q are in reasonable agreement with those obtained 



from STM measurements. p5[ That 7 > 1 for the Gaussian model expresses the fact that 
the He-adsorbate interaction range extends beyond the adsorbate unit-cell. In any inversion 
problem, the question of stability is major importance. In the present case, convergence to 
the fitted values is overall better for the Gaussian model, which also has a smaller deviation: 
the Gaussian model converges to the same values for initial guesses of < O < 5A, whereas 
convergence is obtained for the Lorentzian model only for the smaller range of < ©0 < 3A. 
The initial value of the parameter 7 has almost no influence on the convergence. When 
compared over a large range of k z values (see figure ||), the two models are seen to essentially 
differ in their prediction of the cross section merely by a shift. The situation is different, 
however, when we compare the sought-after quantity: the shape function. The corresponding 
shape functions are shown in figure |9|. The Lorentzian model is much more peaked and 
narrow. Intuitively, it seems that the Gaussian shape is better suited to describe the adatom 
shape function, in agreement with the relative deviations and convergences. 
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D. Inversion of Potential for Weakly Corrugated Surfaces 



Consider a surface of substrate type A, in which one of the atoms is replaced by an atom 
of type B, such that the radii differ only slightly. This is an example of a surface with weak 
corrugation due to substitutional disorder. Let the surface be flat and located at z — 0. The 
interaction potential with an incident He atom may be written approximately as: 

U(R, z) = U.(z) + U d (R 2 + z 2 ) , \U 8 \ > \U d \ (51) 

where R = (x,y), V s = %-U s is the contribution of the surface, and Vd = is the 

interaction with the defect, assumed spherically symmetric. We next assume V s (z) is known 
and wish to determine Vd from the scattering data. We note that a similar inversion has 
already been considered by Gerber, Yinnon and coworkers, |31~1 , [32| but for defectless, crys- 
talline surfaces. When the condition \U S \ ^> \Ud\ is used in eq.(||), a first order expansion 



yields: 



T]0(R) 



poo 

/ dz 


(k 2 -U s {z)f -k z 


Izt(R) 





V (R)^r, (R) + ^ Vl (R) 



- k z z t (R) 



mW—r * Ud(R2 + z2 \ (52) 

J »W (k?-U,(z))> 
The turning-point function Zt(R) may either be measured by the inversion procedure de- 



scribed in section |IV|IVB| for the shape function £, or approximated as the solution to 



V a (z) = E, i.e., neglecting Vd altogether. In either case r] (R) is a known quantity. Let us 



next use our approximation for r\ in eq. ( 34 ) . Then after a first order expansion about 2rj t 



XY- 



ak 2 z Z(k z ,6 in ) --R 2 S + f 3 dRR cos(2r/ ) w / ' ' dRR sm(2r] (R)) r] X {R) (53) 

Zj J J 

The LHS consists of measurable/known quantities and we denote it M\{k z ). The RHS 
contains the unknown function Ud which we are after; we denote it I(k z ). Using the definition 
of 77a (eq.pi)): 



21 



r R s r°° UAR 2 + z 2 ) 

I(k z ) = - dRRsm{2 Vo {R)) dz dK \ (54) 

7o J «m (k 2 -ujz)Y 



A practical approach at this point would be to expand Ud(R 2 + z 2 ) in some convenient basis 
{f n (R 2 + z 2 )}: U d (R 2 + z 2 ) = Ec n fn(R 2 + z 2 ), so that: 

^ r R * r°° f (R 2 + z 2 ) 

Yc n dRR sin (2r] (R) ) / dz JnK ' l = -MAk z ) (55) 

Jo ' J «*> (k z 2 -U s (z)Y 

This linear system can then be solved for as many c n as there are data points (in k z ). A 
convenient choice of (non-orthogonal) basis functions are Gaussians, since then the double 
integral decouples into a product (neglecting the i?-dependence of z t ). This approach should 
give a reasonable approximation to the defect potential. Powerful methods to deal with 
related numerical inversion problems, employing functional sensitivity analysis, have been 



developed by Ho and Rabitz. [35 



However, with some further assumptions it is also possible to perform an actual inversion 
for Ud, by combining the technique of section [V LV B and an Abel transform. We proceed 
to show this. 

Our first two additional assumptions are: 

• The interaction with the surface is purely repulsive, so that ( = U s (z) is a 
monotonously decreasing function. 

• Also p = rjo(R) is a monotone function. That this is reasonable can be seen by 
considering a hard-wall system; then a monotonously decreasing shape function (see 
section |TV[IV B| ) induces a monotonously increasing phase shift, by eq. fl29|) . 



Under these assumptions we may consider the inverse functions z(() and R(p). The 
transformation to the variable ( decouples the integration limits if we make the fair approx- 
imation that k 2 = U a (zt(R)) (i.e., neglecting Ud), for then: 

dz grfgif^/y* wo 



zt[R) {k 2 -u s {z)y h > dC (k 2 -c 
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and hence: 

I(k z ) = £ d( (k z 2 - C)"V(C) where 
dz r^s 

/(C) = — / di2 i2 sin(2r ?0 (i?)) E^fl, C) (56) 
aQ Jo 

The above expression has the form of the familiar Abel transform, for which an explicit 
inversion exists: 

™ (<-*)-' < 57 » 

We rewrite this, using / = Mi (eq.(|53D) and changing variables from R to p, as: 

rP( R °) dR 
J = - / dp C/ d (p, C) i?(p) sin(2p) = G((, k z ) (58) 
Jp(o) dp 

where: 



The function G((, k z ) consists entirely of known quantities. We are still left with the task 
of isolating the potential But it is clear that J in eq . (|58"D is almost identical to Iq in 



eq . (|39D . In order to be able to use the Fourier transform technique applied to Iq, we must 



now introduce one additional assumption: 

• T]o(R) may reasonably well be approximated in the "hard- wall" form: t]q(R) = 
—k z ^(R). Here £(R) is to be considered known from the inversion of cross section 
data as described in section |IVPVB| Hamburger et al. [0] have shown that this form 



holds quite well also for realistic potentials. 
Changing variables to 6 = 2£(R): 

dR 

W.O^^e) sin(^e) 

2 : JqT' d<S> F ^ ( e * fcze ~ e ~ lkz@ ) = Ie ( 60 ) 



j= de 

Je(o) 

1 fO(Rs) 



where now: 
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Repeating the arguments leading to eq.(f45|) and remembering that now I = G((,k z ), we 
finally obtain: 

1 4? 

W> = lavm ^r^ 1 h); 8] (62) 
Formally, this completes the inversion. 



V. CONCLUDING REMARKS 

In this paper we investigated theoretically the cross section of a defect on a surface. 
We found that this quantity, when available from experiments, contains important infor- 
mation on the defect and its interaction with an incoming atom. We demonstrated that 
the cross section can be used for such purposes as an inversion to yield the shape function 
of a defect, an approximate inversion to give the interaction potential with a substituted 
atom, and evaluation of the role of the long-range repulsive forces in the interaction with 
an incident atom. The single most important advantage of the cross section is probably 
that it is so easily measurable, by employing the definition of the cross section proposed 
by Poelsema and Comsa. |l[ This operational definition was shown here to be equivalent 
to the standard formal definition of quantum scattering theory, and hence either can be 
used interchangeably. Using the formal definition and the Sudden approximation, |23] we 



derived a form of the Optical theorem which is quite amenable to analytical study. Future 
work may benefit from this expression for the cross section. For example, it seems that 
extracting various probability distributions characterizing dilute surface disorder (such as 
a distribution of radii for hemispherical defects) should prove possible within this frame- 
work. Traditional approaches tended to rely on differential cross section measurements for 
inversion applications. Our general message for future work is that the total cross section, 
which is much easier to obtain experimentally, will be a highly fruitful subject of study for 
the purpose of both structural and dynamical characterization of isolated defects. Indeed, 
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a very promising possibility for future work should be one of measuring cross sections for 
surface defects of interest, and using this data in a corresponding theoretical effort to extract 
information on the position of the defect, its geometric shape, and its interaction potential 
with the incident atom. The latter is especially important, since the interaction potential 
contains information on the electronic density structure of the defect. We emphasize that 
cross sections measurements, when performed over a wide range of energies and incidence 
angles, may well become a new type of surface microscopy. We are currently pursuing efforts 
along this line, in cooperation with the experiments of P. Zeppenfeld. 
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FIGURES 

FIG. 1. Cross section in the cone case (n = 1), for two values of the corrugation parameter h. 
The cross sections are almost identical, except for the behavior at small k z . 

FIG. 2. Cross section in the convex case (n = 2), for two values of the corrugation parameter h. 
The stronger the corrugation, the faster and more pronounced are the oscillations, and the larger 
the cross section for small k z . 

FIG. 3. Cross section in the convex case, for three values of the convexity parameter n, at a 
constant corrugation h = 1A. The higher the convexity, the faster and more pronounced are the 
oscillations, and the larger the cross section for small k z . 

FIG. 4. Cross section in the concave case for two values of the concavity parameter n, at a 
constant corrugation h = 1A. The larger the concavity, the smaller the cross section. Oscillations 
are virtually unnoticeable. 

FIG. 5. Cross section for k z — > 0, in the concave and convex cases, as a function of the 
concavity /convexity parameter. The more concave (i.e. peaked), the smaller the cross section. 
The opposite happens the more convex (i.e. rectangular) the shape function. 

FIG. 6. General form of the shape function Q(R) and its inverse R(Q). The shape function 
should have a smooth top and monotonously decrease over a substantial range. We ignore the 
possibility of a subsequent increase. 

FIG. 7. Cross section data for single Ag atom on Pt(lll) (circles), along with fits for Lorentzian 
(solid line) and Gaussian (dashed line) shape function models. The agreement with the Gaussian 
model is somewhat better. 

FIG. 8. Behavior of the cross sections of the Lorentzian and Gaussian models over a large range 
of incidence wave numbers. The two models predict a similar behavior, differing essentially by a 
shift along the y-axis. 
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FIG. 9. Shape functions of the Ag adatom for the Lorentzian and Gaussian models, using the 
parameters obtained by fitting the predicted cross sections to the measured one. The Gaussian 
model has a smoother top and is less narrowly peaked, in agreement with intuitive expectation for 
the "correct" shape. Note the different scales along the x and y-axes. 
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